Dirichlet Boundary Conditions Research Articles

Global weak solvability in a self-consistent chemotaxis-Navier-Stokes system involving Dirichlet boundary conditions for the signal

Global weak solvability in a self-consistent chemotaxis-Navier-Stokes system involving Dirichlet boundary conditions for the signal

A Physically Motivated Heat Source Model for Laser Beam Welding

In this contribution, we present a physically motivated heat source model for the numerical modeling of laser beam welding processes. Since the calibration of existing heat source models, such as the conic or Goldak model, is difficult, the representation of the heat source using so-called Lamé curves has been established, relying on prior Computational Fluid Dynamics (CFD) simulations. Lamé curves, which describe the melting isotherm, are used in a subsequent finite-element (FE) simulation to define a moving Dirichlet boundary condition, which prescribes a constant temperature in the melt pool. As an alternative to this approach, we developed a physically motivated heat source model, which prescribes the heat input as a body load directly. The new model also relies on prior CFD simulations to identify the melting isotherm. We demonstrate numerical results of the new heat source model on boundary-value problems from the field of laser beam welding and compare it with the prior CFD simulation and the results of the Lamé curve model and experimental data.

The existence of positive solutions for a critically coupled Schrödinger system in a ball of [formula omitted]

In this paper, we consider the following coupled Schrödinger system with doubly critical exponents, which can be seen as a counterpart of the Brezis-Nirenberg problem{−Δu+λ1u=μ1u5+βu2v3,x∈Ω,−Δv+λ2v=μ2v5+βv2u3,x∈Ω,u=v=0,x∈∂Ω, where Ω is a ball in R3, −λ1(Ω)<λ1,λ2<−14λ1(Ω), μ1,μ2>0 and β>0. Here λ1(Ω) is the first eigenvalue of −Δ with Dirichlet boundary condition in Ω. We show that the problem has at least one positive solution for all 0<β<2μ1μ2. In particular, when λ1=λ2, we prove the existence of positive synchronized solutions for all β>0.

Multinode Shepard method for two-dimensional elliptic boundary problems on different shaped domains

In this paper, we continue the study on the application of multinode Shepard method to numerically solve elliptic Partial Differential Equations (PDEs) equipped with various conditions at the boundary of domains of different shapes. In particular, for the first time, the multinode Shepard method is proposed to solve elliptic PDEs with Dirichlet and/or Neumann boundary conditions. The method has been opportunely handled to efficiently work dealing with scattered distribution of points and, to this aim, several experiments in different 2d domains have been performed. Comparisons with the analytic solution and the results generated by the Kansa’s RBF solvers have been reported referring to Halton points. The results are very promising and should be of interest for applications in the real world.

Efficient atomistic/continuum coupling using lattice Green’s functions

Coupled Atomistic/Continuum (A/C) simulation methods significantly reduce the number of degrees of freedom in the solution of a boundary value problem by using continuum elasticity in domains where defects interact elastically. A/C coupling involves two boundaries: (a) an inner boundary that connects the atomistic and continuum domain and (b) an outer boundary where displacement or force boundary conditions are applied. Here, the lattice Green’s function (LGF)–based flexible boundary condition method is combined with the LGF-based discrete boundary element method to satisfy coupling and boundary conditions at the inner and outer boundaries, respectively, of 3D domains. The coupling algorithms for Dirichlet boundary conditions are presented, with other cases being similar. The LGF is further replaced by a more efficient coarsened LGF method that coarsens the atomistic degrees of freedom at the outer boundary. The scaling of memory and computation versus 3D size and degree of coarsening are discussed. The method and its accuracy are demonstrated using a large 2D problem involving applied forces inside the atomistic domain.

A novel Touchard polynomial-based spectral matrix collocation method for solving the Lotka-Volterra competition system with diffusion

This paper presents the computational solutions of a time-dependent nonlinear system of partial differential equations (PDEs) known as the Lotka-Volterra competition system with diffusion. We propose a combined semi-discretized spectral matrix collocation algorithm to solve this system of PDEs. The first part of the algorithm deals with the time-marching procedure, which is performed using the well-known Taylor series formula. The resulting linear systems of ordinary differential equations (ODEs) are then solved using the spectral matrix collocation technique based on the novel Touchard family of polynomials. We discuss and establish the error analysis and convergence of the proposed method. Additionally, we examine the stability analysis and the equilibrium points of the model to determine the stability condition for the system. We perform numerical simulations using diverse model parameters and with different Dirichlet and Neumann boundary conditions to demonstrate the utility and applicability of our combined Taylor-Touchard spectral collocation algorithm.

Stability of steady‐state solutions of a class of Keller–Segel models with mixed boundary conditions

In this paper, we investigate the existence and stability of non‐trivial steady‐state solutions of a class of chemotaxis models with zero‐flux boundary conditions and Dirichlet boundary conditions on a one‐dimensional bounded interval. By using upper–lower solution and the monotone iteration scheme method, we get the existence of the steady‐state solution of the chemotaxis model. Moreover, by adopting the “inverse derivative” technique and the weighted energy method, we prove the stability of the steady‐state solution of this chemotaxis model.

Renormalized stress-energy tensor on global anti-de Sitter space-time with Robin boundary conditions

We study the renormalized stress-energy tensor (RSET) for a massless, conformally coupled scalar field on global anti-de Sitter space-time in four dimensions. Robin (mixed) boundary conditions are applied to the scalar field. We compute both the vacuum and thermal expectation values of the RSET. The vacuum RSET is a multiple of the space-time metric when either Dirichlet or Neumann boundary conditions are applied. Imposing Robin boundary conditions breaks the maximal symmetry of the vacuum state and results in an RSET whose components with mixed indices have their maximum (or maximum magnitude) at the space-time origin. The value of this maximum depends on the boundary conditions. We find similar behaviour for thermal states. As the temperature decreases, thermal expectation values of the RSET approach those for vacuum states and their values depend strongly on the boundary conditions. As the temperature increases, the values of the RSET components tend to profiles which are the same for all boundary conditions. We also find, for both vacuum and thermal states, that the RSET on the space-time boundary is independent of the boundary conditions and determined entirely by the trace anomaly.

Optimal regularity and fine asymptotics for the porous medium equation in bounded domains

Abstract We prove the optimal global regularity of nonnegative solutions to the porous medium equation in smooth bounded domains with the zero Dirichlet boundary condition after certain waiting time T * {T^{*}} . More precisely, we show that solutions are C 2 , α ⁢ ( Ω ¯ ) {C^{2,\alpha}(\overline{\Omega})} in space, with α = 1 m {\alpha=\frac{1}{m}} , and C ∞ {C^{\infty}} in time (uniformly in x ∈ Ω ¯ {x\in\overline{\Omega}} ), for t &gt; T * {t&gt;T^{*}} . Furthermore, this allows us to refine the asymptotics of solutions for large times, improving the best known results so far in two ways: we establish a faster rate of convergence O ⁢ ( t - 1 - γ ) {O(t^{-1-\gamma})} , and we prove that the convergence holds in the C 1 , α ⁢ ( Ω ¯ ) {C^{1,\alpha}(\overline{\Omega})} topology.

On a class of doubly nonlinear evolution equations in Musielak–Orlicz spaces

AbstractThis paper is concerned with a parabolic evolution equation of the form , settled in a smooth bounded domain of , , and complemented with the initial conditions and with (for simplicity) homogeneous Dirichlet boundary conditions. Here, stands for a diffusion operator, possibly nonlinear, which may range in a very wide class, including the Laplacian, the ‐Laplacian for suitable , the “variable‐exponent” ‐Laplacian, or even some fractional order operators. The operator is assumed to be in the form with being measurable in and maximal monotone in . The main results are devoted to proving existence of weak solutions for a wide class of functions that extends the setting considered in previous results related to the variable exponent case where . To this end, a theory of subdifferential operators will be established in Musielak–Orlicz spaces satisfying structure conditions of the so‐called ‐type, and a framework for approximating maximal monotone operators acting in that class of spaces will also be developed. Such a theory is then applied to provide an existence result for a specific equation, but it may have an independent interest in itself. Finally, the existence result is illustrated by presenting a number of specific equations (and, correspondingly, of operators , ) to which the result can be applied.

Casimir and Helmholtz forces in one-dimensional Ising model with Dirichlet (free) boundary conditions

Attention in the literature has increasingly turned to the issue of the dependence on ensemble and boundary conditions of fluctuation-induced forces. We have recently investigated this problem in the one-dimensional Ising model with periodic and antiperiodic boundary conditions (Annals of Physics 459, 169533 (2023)). Significant variations of the behavior of Casimir and Helmholtz forces was observed, depending on both ensemble and boundary conditions. Here we extend our study by considering the problem in the important case of Dirichlet (also termed free, or missing neighbors) boundary conditions. The advantage of the mathematical formulation of the problem in terms of Chebyshev polynomials is demonstrated and, in this approach, expressions for the partition functions in the canonical and the grand canonical ensembles are presented. We prove analytically that the Casimir force is attractive for all values of temperature and external ordering field, while the Helmholtz force can be both attractive and repulsive.

Simulation of viscoelastic free-surface flows with the Particle Finite Element Method

AbstractViscoelastic fluids are central in numerous applications from polymer manufacturing to the pharmaceutical industry and biological research. However, since analytical solutions are generally not available or too complex, it is common practice to study free-surface viscoelastic flows through numerical simulation techniques. This work proposes the use of the so-called particle finite element method (PFEM), a Lagrangian approach combining standard FEM techniques with a remeshing strategy. The PFEM is able to efficiently handle mesh distortion and to accurately track the free-surface evolution. Therefore, it is exploited in this work to deal with large displacements problems in the context of nonlinear viscoelasticity. An implementation of the Oldroyd-B constitutive model in the PFEM framework is here presented including details regarding how to deal with the transfer of the internal variables during remeshing events. Additionally, an innovative approach to impose unilateral Dirichlet boundary conditions ensuring optimal mass conservation is presented. The implementation is verified with two free-surface highly viscous benchmark flows: the impacting drop and the jet buckling problems. The results show perfect agreement with those obtained with other numerical techniques. The proposed framework opens the way for using PFEM in various applications, ranging from polymer extrusion to more sophisticated scenarios involving viscoelastic and viscoelasto-plastic constitutive laws.

Exponential decay of solutions of damped wave equations in one dimensional space in  the $L^p$ framework for various boundary conditions

We establish the decay of the solutions of the damped wave equations in one dimensional space for the Dirichlet, Neumann, and dynamic boundary conditions where the damping coefficient is a function of space and time. The analysis is based on the study of the corresponding hyperbolic systems associated with the Riemann invariants. The key ingredient in the study of these systems is the use of the internal dissipation energy  to estimate  the difference of solutions with their mean values in an average sense.

Qualitative analysis on a reaction–diffusion SIS epidemic model with nonlinear incidence and Dirichlet boundary

In this paper, the dynamical behavior in a spatially heterogeneous reaction–diffusion SIS epidemic model with general nonlinear incidence and Dirichlet boundary condition is investigated. The well-posedness of solutions, including the global existence, nonnegativity, ultimate boundedness, as well as the existence of compact global attractor, are first established, then the basic reproduction number R0 is calculated by defining the next generation operator. Secondly, the threshold dynamics of the model with respect to R0 are studied. That is, when R0<1 the disease-free steady state is globally asymptotically stable, and when R0>1 the model is uniformly persistent and admits one positive steady state, and under some additional conditions the uniqueness of positive steady state is obtained. Furthermore, some interesting properties of R0 are established, including the calculating formula of R0, the asymptotic profiles of R0 with respect to diffusion rate dI, and the monotonicity of R0 with diffusion rate dI and domain Ω. In addition, the bang–bang-type configuration optimization of R0 also is obtained. This rare result in diffusive equation reveals that we can control disease diffusion at least at one peak. Finally, the numerical examples and simulations are carried out to illustrate the rationality of open problems proposed in this paper, and explore the influence of spatial heterogeneous environment on the disease spread and make a comparison on dynamics between Dirichlet boundary condition and Neumann boundary condition.

Asymptotic Errors in the Superconvergence of Discontinuous Galerkin Methods Based on Upwind-Biased Fluxes for 1D Linear Hyperbolic Equations

AbstractIn this paper, we study the superconvergence of the semi-discrete discontinuous Galerkin (DG) method for linear hyperbolic equations in one spatial dimension. The asymptotic errors in cell averages, downwind point values, and the postprocessed solution are derived for the initial discretization by Gaussian projection (for periodic boundary condition) or Cao projection Cao et al. (SIAM J. Numer. Anal. 5, 2555–2573 (2014)) (for Dirichlet boundary condition). We proved that the error constant in the superconvergence of order $$2k+1$$ 2 k + 1 for DG methods based on upwind-biased fluxes depends on the parity of the order k. The asymptotic errors are demonstrated by various numerical experiments for scalar and vector hyperbolic equations.

Assessment of Local Radial Basis Function Collocation Method for Diffusion Problems Structured with Multiquadrics and Polyharmonic Splines

This paper aims to systematically assess the local radial basis function collocation method, structured with multiquadrics (MQs) and polyharmonic splines (PHSs), for solving steady and transient diffusion problems. The boundary value test involves a rectangle with Dirichlet, Neuman, and Robin boundary conditions, and the initial value test is associated with the Dirichlet jump problem on a square. The spectra of the free parameters of the method, i.e., node density, timestep, shape parameter, etc., are analyzed in terms of the average error. It is found that the use of MQs is less stable compared to PHSs for irregular node arrangements. For MQs, the most suitable shape parameter is determined for multiple cases. The relationship of the shape parameter with the total number of nodes, average error, node scattering factor, and the number of nodes in the local subdomain is also provided. For regular node arrangements, MQs produce slightly more accurate results, while for irregular node arrangements, PHSs provide higher accuracy than MQs. PHSs are recommended for use in diffusion problems that require irregular node spacing.

Fourier series method for the stability solution of simply supported thin beams on two-parameter elastic foundations of the Pasternak, Filonenko-Borodich, Hetenyi or Vlasov models

The analysis of stability problems of beams on two-parameter foundations (Bo2PFs) is an important part of their design for compressive loads. This work presents novel first principles derivation of the governing differential equations of elastic stability (GDES) of thin beams resting on two-parameter elastic foundations of the Pasternak, Filonenko-Borodich, Hetenyi or Vlasov models. The requirements of translational and rotational equilibrium of all the applied, reactive and internal forces on an infinitesimal segment of the Bo2PF and the laws of infinitesimal calculus were used to formulate the GDES as a fourth order ordinary differential equation (ODE) in terms of the transverse displacement function ux. The GDES is non-homogeneous in the presence of applied transverse load qx but homogeneous when qx vanishes. This study presents the Fourier series method (FSM) for solving the governing differential equation of stability (GDES) for the case of Dirichlet boundary conditions. The FSM has the advantage of amenability to differentiation, and integration due to the orthogonality properties of the sinusoidal functions. Implementation of the FSM by assuming the unknown function in the GDES as a Fourier series of infinite terms and the exploitation of orthogonalization simplifies the problem to an algebraic eigenvalue problem which is the characteristic buckling equation. The exact eigenvalues are found by algebraic solution of the buckling equation. The exact eigenvalues were used to find the exact buckling loads and the exact buckling load coefficients. The critical buckling load was found to correspond to the first buckling mode (n= 1), and is identical with previous solutions in the literature. Numerical calculations for the critical buckling load parameters Kcr were presented for the Bo2PF problem for values of the dimensionless foundation parameters k-1= 0, k-2= 0; k-1= 100, k-2= 0; k-1= 0, k-2= 1; k-1= 100, k-2= 100; k-1= 0, k-2= 2.5; k-1= 100, k-2= 2.5. The present solutions were compared with previous solutions for Kcr in the literature. The comparison shows that the present FSM results are identical with previous results obtained using various other methods such as Recursive Differentiation Method, Finite Element Method, Generalized Integral Transform Method (GITM) and Stodola-Vianello Iteration Method. The study has illustrated the effectiveness of the FSM for solving Bo2PFs.

On the space-time analyticity of the inhomogeneous heat equation on the half space with Neumann boundary conditions

We consider the inhomogeneous heat equation on the half-space R + d with Neumann boundary conditions. We prove a space-time Gevrey regularity of the solution, with a radius of analyticity uniform up to the boundary of the half-space. We also address the case of homogeneous Robin boundary conditions. Our results generalize the case of homogeneous Dirichlet boundary conditions established by Kukavica and Vicol [A direct approach to Gevrey regularity on the half-space. In: Partial differential equations in fluid mechanics. Cambridge: Cambridge Univ. Press; 2018. p. 268–288. (London math. soc. lecture note ser.; vol. 452).].

Strong solutions and attractor dimension for 2D NSE with dynamic boundary conditions

AbstractWe consider incompressible Navier–Stokes equations in a bounded 2D domain, complete with the so-called dynamic slip boundary conditions. Assuming that the data are regular, we show that weak solutions are strong. As an application, we provide an explicit upper bound of the fractal dimension of the global attractor in terms of the physical parameters. These estimates comply with analogous results in the case of Dirichlet boundary condition.

Analysis of the boundary conditions for the ultraweak-local discontinuous Galerkin method of time-dependent linear fourth-order problems

In this paper, we study the ultraweak-local discontinuous Galerkin (UWLDG) method for time-dependent linear fourth-order problems with four types of boundary conditions. In one dimension and two dimensions, stability and optimal error estimates of order k + 1 k+1 are derived for the UWLDG scheme with polynomials of degree at most k k ( k ≥ 1 k\ge 1 ) for solving initial-boundary value problems. The main difficulties are the design of suitable penalty terms at the boundary for numerical fluxes and the construction of projections. More precisely, in two dimensions with the Dirichlet boundary condition, an elaborate projection of the exact boundary condition is proposed as the boundary flux, which, in combination with some proper penalty terms, leads to the stability and optimal error estimates. For other three types of boundary conditions, optimal error estimates can also be proved for fluxes without any penalty terms when special projections are designed to match different boundary conditions. Numerical experiments are presented to confirm the sharpness of theoretical results.